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1^ Abstract 

O 

We develop the holographic formulation of transport in strongly coupled two- 
layer systems. We identify a dc conductivity, a dc , that is finite even in a 

> 

translationally invariant setup, and universal for CFTs with a gravity dual. 

The thermoelectric conductivity and heat conductivity are fully determined 

by the electrical conductivity matrix, as a consequence of Ward identities. We 

use the memory-matrix approach for double-layer systems, together with Ward 

identities, to show that a dc - extended to finite frequency - has no Drude peak 

and, similarly, that its universal value is unaffected if translation invariance is 

softly broken. 



1 Introduction 

Charge transport is a fundamental manifestation of the underlying physical principles of 
any condensed-matter system and thus a quantity of intense interest in any given theory 
or experiment. In a field-theory approach, transport is characterized by near-equilibrium 
correlation functions with lorentzian signature. At strong coupling these correlation func- 
tions are often inaccessible, due to the reliance on perturbative methods and due to severe 
restrictions of lattice simulations for finite-density quantities that are intrinsically defined 
in lorentzian signature. It is generally accepted as an attractive feature of the AdS/CFT 
correspondence [U El E] that it provides a framework for the computation of transport 
quantities for strongly- interacting field theories [4J. This approach is most promising for 
transport coefficients that are universal within models with a gravity dual. The most well- 
known such transport coefficient is the ratio of shear viscosity over entropy density, rj/s 
of a relativistic strongly coupled plasma with a gravity dual, which takes on the universal 
value of 1/4-7T in suitable units [3j. Such quantities, however, are hard to come by. It is of 
great interest and importance to continue the search for transport coefficients and univer- 
sal mechanisms if we are to push the applicability of holographic methods further towards 
realistic systems. Applying the techniques of holography to condensed matter phenom- 
ena is a promising arena in which to look for such universal transport effects. There is, 
however an immediate obstruction when we consider the properties of theories at finite 
charge density. The most directly accessible quantities, namely dc charge-transport co- 
efficients, are ill defined in the continuum limit, since they are strictly divergent. The 
divergence is directly linked to translation invariance; its breaking - either by impurities 
or via lattice effects - is likely not universal. It is nevertheless sometimes possible, for 
example by exploiting the hydrodynamic nature O El E] of the late-time transport, to 
extract scaling laws [9j [TOj [Til H2] for dc conductivities. These mechanisms rely on a finite 
Drude weight in the conductivity, associated with approximately conserved quantities. In 
this paper we take a different angle, focusing on a quantity with vanishing Drude weight 
and derive universal results which do not involve hydrodynamic modes and can therefore 
be considered intrinsic to the microscopic field theory. 

To this end we introduce a new transport coefficient, trans conductance - or its inverse 
transresistance - into the holographic dictionary and show that a suitably defined transport 
coefficient characterizing transconductance is finite even in the translationally invariant 
setup, and in fact turns out to take a universal value for (2+1) dimensional CFTs. The 
dc transconductance exhibits complete bulk universality and has the value 

a dc = -^, (1.1) 

e eff 



in any 2+1 dimensional quantum field theory with a gravity dual that asymptotes to AdS±. 
The effective bulk coupling e^ s is related to the number of charged degrees of freedom, 
i.e. a central charge, of the IR fixed point of this theory. There is a well-known exception 
to our earlier statement that dc charge transport must be infinite in a translationally 
invariant setup [13]: a particle-hole symmetric system (such as a zero-density CFT) can 
give rise to a finite current, because a finite charge current involving particles and holes 
in equal numbers carries no net momentum. The transport we are describing here is a 
realization of an analogous mechanism, but applies away from the particle-hole symmetric 



point. In other words the transport associated with the universal coefficient (1.1 ) involves 
a charge flow with no associated net momentum flow and is thus finite for a clean system. 
As far as we know this is the first time such a mechanism is pointed out in a system that 
has a net charge density and is not at the particle-hole symmetric point. It likely has 
implications for other systems and approaches, holographic or not. 

One can furthermore use the memory matrix approach to show that the universal value 
of <r dc only receives parametrically small corrections if translation invariance is broken via 
a coupling to the leading irrelevant symmetry-breaking operator, ensuring that the theory 
still flows to the same IR fixed point. Although our approach applies to quantum critical 
dynamics, we expect our results to be more widely valid, indeed they may well follow from 
an analysis of the constraints placed upon our system by relativistic hydrodynamics, for 
example within a Boltzmann approach [13]. It has been suggested that clean graphene 
monolayers can be thought of as a relativistic critical system [H] [15], at least for a 
window of temperatures where the effective speed of light runs very little with scale, so 
our approach to the physics of double layers may be of direct relevance to the case of 
graphene. For a probe-brane holographic approach to double- layer graphene see [T6] . 

This paper is structured as follows. In the remainder of this introduction we give a def- 
inition of transconductance, the transport at the focus of this work. Section [2] introduces 
the general class of holographic bottom-up models which exhibit transconductance and 



uses general arguments to establish Eq. (1.1). Section 3 serves to anchor a in a more 



general discussion of the conductivity matrix of the models introduced in Section [2] and 
then discusses constraints placed on transport by Ward identities. Section [4] considers 
the effect of translation breaking via impurities or lattice effects and uses the memory 
function matrix approach (with a short review) together with the clean Ward identities 
to establish general forms of low-frequency limits of transconductance. The conclusions 
are followed by a short appendix on some details of the holographic renormalization of 
our action. 



1.1 Introducing: a new transport coefficient 

Transresistance, conceptually, measures the amount of drag between charge carriers in two 
different layers^] under the condition that there is no direct interaction between them. 
Systems where such a quantity can be investigated typically have two conducting layers 
separated by an insulating layer, for example two sheets of graphene sandwiching a sheet 
of boron nitride [17J. One often refers to the whole setup as a hetero structure. The case 
where the spacing between the layers is so small that quantum-mechanical wave-function 
spread means that the whole heterostructure can be viewed as a two-dimensional system 
is particularly interesting and is the focus of this work. 

Drag is defined as follows: One measures the induced voltage in one layer due to an 
applied current in the other. More precisely, one applies a current J" in layer one (the 
'active' layer), and then measures the induced voltage V% in layer two, the 'passive' layer, 
under the condition that there is no current flowing in the passive layer. This phenomenon 
is often referred to as Coulomb drag, since in many systems, such as graphene heterostruc- 
tures or Hall double layers, interlayer interactions are mediated by the Coulomb force (see 
Fig. §. 

In a field-theory treatment of this phenomenon the basic object which one needs to 
determine is the conductivity matrix a^ b , defined in terms of a Kubo formula [TBI US] as 



Lib 



:«,Q) = -^ydte(t)e tot ([jJ(*,Q) ) JJ(0,Q)]>. (1.2) 



In this expression A is the area of the sample, i,j are spatial indices and a, b label the 
individual layers. The integrand is nothing other than the retarded correlation function 
of the conserved currents in both layers 

G% R (t,Q) = -iO(t){[Ji(t,Q),JmQ)])- (1-3) 

In this paper we use holography to calculate this correlation function and thus the 
transconductance in a strongly-coupled double-layer system. The conductivity associated 
with a single U(l) sector was the object of much attention (see e.g. [201 EB [221 [231 121]) 
and we shall use these results as a benchmark for the ones here. In this paper we are 
interested in the conductivities themselves and so we will mostly study the correlation 



functions at zero momentum Q = 0. Of course the full correlation function (1.2) contains 
more physical information, which will be useful, for example, in developing a hydrody- 
namic description of drag. In perturbative treatments of drag, one often has the restriction 



■"■We are interested in this paper in two-dimensional sheets or layers, each governed by a 2+1 strongly- 
interacting field theory. 
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Figure 1: a) The leading contribution to drag in a holographic computation. 
The two currents interact via mixing with the momentum operator Tq x := P. 
The diagram contributing to the finite quantity a involves no net momentum 
transfer, as indicated by the two opposite momentum-flow arrows, b) Example 
of contribution in perturbative formulation of drag. The currents have no direct 
interaction, but talk to each other via loops in each layer that in turn couple through 
e.g. the Coulomb interaction (wiggly lines). The diagram with only a single wiggly 
line joining the two fermion loops vanishes at Q = 0. 



012 <C en, <722 ( see e -g- [19J ) , i.e. one assumes that the interlayer drag is very weak. In 
this work we shall have no such restriction, and we are able to analyze the case where 
the drag itself is as strong as the intralayer interactions. These interactions are those of 
a strongly-coupled field theory. 



2 Finite dc transport with translation invariance 

From the discussion above we can see that the requirements on our model are simple. 
Each of the conserved currents in the theory gives rise to a local U(l) gauge field Am 
in the bulk. Since we do not want to consider situations where the two currents have a 
direct interaction, we do not include a cross kinetic term for the two gauge fields in the 
bulk, which in principle would be allowed by the U(l) x U(l) local symmetry. Generic 
supergravity models usually have gauge kinetic terms parametrized by a (neutral) scalar 
field - we refer to this as a 'dilaton' even though it is not the stringy dilaton field. Hence 
we also include a field <p with a potential V(<p). We shall not need the exact form of the 
potential or the function $[y>] that governs the coupling of the dilaton to the gauge fields, 
for our argument, but we should require that the field p has at least marginal fall off as 
one approaches the boundary and that, $(z) — ¥ 1 in this limit. Indeed in this paper, we 
are concerned with generic quantities and mechanisms, which will be necessarily satisfied 
in any specific microscopic realization. Since we consider both of our layers as equivalent 
physical systems individually, their conserved currents necessarily couple equally to the 



rest of the theory. As a consequence the effective bulk description is of the form 

s = / ^ (h { R + 1 ~ m2 ~ 2V ^) ~ is ^> + F «>) d ' x (21) 

and we see that it has an SO (2) symmetry rotating the gauge fields into one another. We 
can think of this as a flavor symmetry of currents. If only a layer-exchange Z2 symmetry 



remains, one could allow more general couplings than those in (2.1), in particular each 
layer could then couple differently to the remaining degrees of freedom. Nevertheless, the 
remaining exchange symmetry is still sufficient to argue finiteness of a dc , as we shall see 



in section 4J3 Holographic theories with two bulk Maxwell fields have been considered 
previously both with 'conventional' AdS asymptotics [251 [2S], arid with a few on non- 
relativistic holography 



In section 3.2 we will add a suitable counterterm action, which makes the variational 



principle well defined and holographically renormalizes the dual theory. The action (2.1) 
does not contain any couplings between Fm and Fr 2 -\, so one might think that there can 
be no interlayer transport. However, as we shall see, interlayer transport is mediated by 
a mixing of the currents with the momentum operator P = T tx , which is holographically 
dual to a fluctuation of the metric, Sg tx . Both gauge fields are coupled to the metric, 
and so the leading-order contribution to drag in our setups proceeds via a momentum 
exchange (graviton exchange) as illustrated in Fig [Ik. Before turning to the computation 
of the conductivities, however, we need to establish some facts about the background 
solutions appropriate for the physics we are interested in. 

2.1 Two-charge solutions and their flavor rotation 

The equilibrium ensemble of the double-layer system is described by a two-charge (dilaton) 



black hole, a solution of the theory (2.1 ) of the form 

<p = <p(z), 

A® = (j)(i){z)dt, z = l,2 

f 2 
ds 2 = G^dsfdx" = — (-f(z)dt 2 + f{z)- l dz 2 + dx 2 ) . (2.2) 

Explicit examples of AdS dilaton black holes [28l 129] , and their IR limits [30] have recently 
been considered in the holographic literature but for the purpose of the analysis here it 
is sufficient to suppose that a solution can be constructed (numerally or analytically). 
Our arguments are made on a formal level and do not rely on any particular solution. 
Consider now the gauge-field equations 

V M ($HF ( O=0. (2.3) 



Plugging in the ansatz (2.2) one deduces that these have a first integral 

C, 



(<&W (0 )' = <j>' {i) {z) 



'(0 



(2.4) 



This can be integrated to give the formal solution 

" z C(i)dz 



(0 



:-h 



§(z) 



(2.5) 



where the lower limit of integration is chosen so that A t (zh) = 0. It is useful to define the 
dimensionless integral 

- 1 dc 



(C = z/z h ) 



(2.6) 



/o <*K0 

which integrates the flow of the coupling constant as a function of the bulk coordinate z. 
This integral may have an infra-red divergence depending on the details of the solution. 
If the dilaton coupling function is trivial, we have simply 1$ = 1. The constants Cu\ are 
the charge densitieqj associated with the gauge fields, as can be seen from the integral 

1 /"_,,_ C (i) 



(rii) 



2e 2 A 



*®[<p]F(i) 



2e 2 



(2.7) 



where * denotes the Hodge dual and A is the area of the spatial IR 2 of the dual field theory 



as in (1.2). This is equivalent to defining the charges by differentiating the action with 



respect to the boundary value of <j)u\. The boundary limit of the expression (2.5) in fact 
defines the chemical potentials, which can be brought into the form 



Hi = 2e UiZhU ■ 

This allows us to determine the charge susceptibility Xij 
computes 

X(T) 



drij 



ii=0 



(2.8) 
x(T)Sij. One 

(2.9) 



2e 2 z h (T)U ' 

The dependence on temperature is left implicit and is hidden in the horizon radius Zh{T), 
evaluated at vanishing charge density. If we now define the quantity p, = C( 2 )/C(i) we can 
construct the SO (2) matrix 



Mifi) 




(2.10) 



2 The sign is convention, coming from the fact that in AdS/CFT the gauge field is dual to a boundary 
current (•//"•>) ~ FL so that the charge density (J? t \) comes with a timelike index up accounting for the 



negative sign in (2.7 1 




which rotates the gauge fields into a new basis 

and acts on the charges as C + = Cm(l + p 2 ) 1 ^ 2 and C_ = 0. The action on the field 
strengths similarly gives 

C 

F_=0, F + = —t-dzAdt. (2.12) 

$(z) 

Hence the black hole is charged under only one combination of the fields in the new basis. 
In order to determine charge transport of this system, we need to solve the quadratic 



fluctuation equations, expanded around the background (2.2). A little thought then 
shows that the fluctuations of A^ cannot couple to any other fluctuations of the system, 
i.e. to fluctuations of the metric, the other gauge field, or the dilaton, for the simple 
reason that these couplings, for example on the level of the action, have to be of the form 
F^5A_5a. We have denoted a generic fluctuation of the other fields as 5a. Since the 
background has F_ = 0, no such couplings can arise. We emphasize that this result is 
completely general, i.e. it holds for a general perturbation Sa(t, z, x) with arbitrary time 
and space dependence. 

2.2 Universal charge transport and its finite dc limit 

In order to determine the dc transport of the model we need to determine the two-point 
correlation functions of the current dual to the gauge fields A^y For this one needs the 
quadratic fluctuation action. We can summarize the discussion above by stating that we 
are now faced with solving the charge transport properties of the current J_ , dual to the 
bulk field A- with an effective theory of the form 

S eS = -^J^^]F 2 + S[g,F + ^}, (2.13) 

which experiences a background that is neutral with respect to the field A_. This implies 
that quadratic fluctuations of A_ completely decouple from the ones of S[g, F + , tp] and we 



only need to consider the first term in (2.13). It then follows directly that the membrane- 



paradigm results of [31] can be applied to the theory at hand. In the present context it 
means that there is a universal dc conductivity associated with A_ , namely 



a dc = — -, (2.14) 

e I {zi R ) 



where one defines the scale-dependent gauge coupling e 2 (z) = &(z)/e 2 and zir is the 
location of the infrared horizon. Alternatively one can derive this result directly from the 



fluctuation equations following from (2.13). What matters is that the answer is universal 
from the perspective of the bulk dual. Unfortunately, it is not model-independent from 
the point of view of the dual field theory due to the dependence on the dilaton profile 
through the scale-dependent coupling e(z). However, for a large class of cases it is in fact 
universal from the field-theory perspective, namely for any CFT3, where it is determined 
by the zero-temperature current- current correlator [3"2"l |2"3"] 



(J,(*),J,(0)) = ^(<V-2^), (2.15) 

with k = -|-2, so that <7 dc (CFT 3 ) = ^. The quantity k is a dimensionless central charge 
of the CFT for charged degrees of freedom. This is a general result, since we did not have 
to use any of the details of the solution of (2.2), which will strongly depend on whatever 
theory we want to solve. Moreover, such a solution, typically will not be analytically 
tractable. Even if the dilaton has a nontrivial profile, the effective coupling constant 
e 2 { z m) can still be related to a central charge at the IR fixed point, but in this case the 
central charge may depend on the chemical potentials and the correlator ( |2.15[ ) will have 
to be evaluated in the IR theory. Dividing the conductivity by the charge susceptibilit)Q 
we obtain the simple expression 

cr dc 

— = 2z h (T)$(z h )I*. (2.16) 

X 

We have thus identified a transport coefficient, a dc which is finite in the dc limit, without 
having to introduce a mechanism of momentum relaxation. Furthermore, we showed that 
the decoupling of the dynamics of the field dual to the transport current giving rise to a dc 
is completely general, and in particular extends to finite-momentum fluctuations. This 
also implies that the value of <r dc is unaffected by the detailed mechanism of momentum 
relaxation, such as weak impurity scattering or lattice effects, simply because the current 
operator J_ has zero overlap with the momentum operator. We shall return to this in 
later sections. 

Finally we should note that the derivation above focused on the dc conductivity, but in 



the case of certain CFTs with a bulk EM duality [2D] , the result (2.14) in fact extends to 
the optical conductivity, meaning that it is independent of frequency and takes the form 

^ drag M = ^. (2.17) 



3 This ratio has previously been obtained by [23 for a theory without dilatonic couplings and our 
result agrees with theirs in this case. 



Relations of this kind have been noted earlier for the optical conductivity of a single 
current, for example by [23J, |20]. In such cases, however, one actually gets a divergent 
dc limit due to translation invariance unless the background is neutral. Then, since the 
current has a non- vanishing overlap with the momentum operator, the mechanism to relax 
momentum will not only resolve this divergence into a Drude peak, but also destroy the 
frequency independence of the optical conductivity. This is in contrast with our case, 
where the relevant current operator has zero overlap with the momentum operator even 
at finite density. Including higher-derivative terms in the bulk action, breaking the self- 
duality of the Maxwell term, is also be expected to modify the frequency-independence 
of the optical conductivity, as was shown in the case of a single current in |33j . 

3 The (rest of the) conductivity matrix 

We saw in the section above that the difference of the fluctuations A_ sees an uncharged 
black hole, while the average of the fluctuations A + couples to the sum of the charges. 
As a consequence of this decoupling the retarded current-currenty correlation functions of 
the dual field theory are diagonal in this basis : 

([J ± ,J ± ]> = Gl ± (3.1) 

We would like to relate this to the conductivity matrix, a^, as defined in the original 
basis. It is easy to convert the corresponding expressions: 

1 1 



Oil 



1 + /S* w 



«[J + ,J + ]> +/ z 2 ([J_,J_]» 



°22 = --^-Ou 2 ([J + ,J + ]> + ([J_,J_]>) . (3.2) 

One can see that each element of the conductance matrix in this basis gets contributions 
both from J + and from J_, and each matrix element thus picks up a dc divergence from 
the J + part. Moreover, our arguments above established that the universal transport 
quantity in this system is associated with the ([«/_, J-}) correlator. This is a physically 
accessible quantity, namely 

a d ^(co) = a 12 (oo)-ft- 1 a ll (uj). (3.3) 



4 Note that there are other non-vanishing correlators, namely those between stress-tensor and current 
([Tox,J±]) — Gt± an d stress-tensor and stress-tensor ([T 0x ,T 0x ]) — G§ T . We will see that they are 
entirely determined in terms of the current-current correlations functions. 



so that 

a dc = lima drag (w). (3.4) 

This definition makes physical sense. Suppose we want to learn about the amount of drag 
between the layers by measuring the 'bare' coefficient a%2- This will get contributions both 
from the momentum transferred between the layers (the drag we are interested in) and 
the momentum dissipated due to interlayer momentum relaxation. If no such interlayer 
momentum relaxation is effective, we instead measure a divergence in the dc transport 
associated with the interlayer translation invariance. It is thus natural to focus on the 



subtracted quantity (3.3). Of course, the entire point of section 2.2 was that there is 



an independent reason to focus on (3.3), namely that it quantifies a universal transport 
property of strongly coupled field theories with gravity duals. Typically one is interested 
in the equal charge case, p, = 1, so that the definition <T drag = cr 12 — er n is independent of 
the densities and Oyi = ^j (G — C++) etc. are particularly simple. 

3.1 Global symmetries and gauge fixing 

In linear response the conductivity is encoded in the retarded correlation function of 



currents at zero momentum via the Kubo formula (1.2) 



k igg> , Q = 0) 
°abM = ~. • (3.5) 

As indicated the correlation function is evaluated at vanishing spatial momentum. In 
situations where parity is broken, as for example in the presence of a magnetic field, there 
can be an off-diagonal Hall contribution to the conductivity. The remaining continuous 
symmetries at Q = then dictate 

c7^ b = a ij S ab + E ij e ab , (3.6) 

where a l i is the direct conductivity and S u is the Hall contribution. Off-diagonal elements 
in layer indices of the conductivity matrix are referred to as trans- or drag- conductivity. 
We have already seen that that <7i2 = o^i- By the transverse SO (2) symmetry we can 
calculate (T12 by considering fluctuations only in the x direction^ At vanishing spatial 
momentum the spin-0 dilaton fluctuation cannot mix to linear order with the spin-1 fluc- 



tuations of the metric which we need to consider. In order to calculate the quantity (3.5) 



we thus only consider vector fluctuations of the metric and gauge fields, of which there 



5 In the absence of a magnetic field, or another mechanism of breaking P, the Hall transconductance 
vanishes. 
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are four. The deviations of the metric and gauge fields from their vanishing background 



(see (2.2)) are parametrized by 

{A^(t,z),g tx (t,z),g zx (t,z)} . (3.7) 

More precisely these functions are related to small fluctuations of metric and gauge field 
as follows: 

I 2 
5G tx (t,z) = —g tx (t,z) 

z A 

I 2 
SG zx (t,z) = —g zx (t,z) 

z A 

8A®(t,z) = AVfrz). (3.8) 

We will work with their Fourier transforms 

Af{t,z) = J^e-^Afi^z) , g tx (t,z) = J ^e^ g tx {u> , z) . (3.9) 

For the problem at hand, it turns out that a convenient choice of gauge sets g zx (u, z) = 0. 
This 'radial' gauge does not fix the diffeomorphism symmetry completely, so that there 
is a residual gauge transformation 

C i9tx = -lut? , C^Af = C igzx = (3.10) 

for constant gauge parameter £ x . This takes the form of an infinitesimal translation 
in the x direction in the field theory and, as we will see later, is tied to translation 
invariance of the boundary theory. Note we only need to concern ourselves with the gauge 
fixing/freedom in the vector sector of which we have just given a complete description. 
Obviously there is more freedom in the spin-0 and spin-2 sectors which do not mix with 
the quantities here at zero momentum. 

3.2 Boundary action, correlation functions & Ward identities 

Correlation functions in AdS/CFT can be calculated by identifying the on-shell action 
as the generating functional of correlation functions in the dual theory. In our case, the 
on-shell action, suitably renormalized, reduces to a boundary term. Before we can identify 
the boundary action we must specify the counter term action further. We take the usual 
form [331 E5] 

Sct = \ I V^7 it + I) d 3 x + Sg T . (3.11) 



k 2 j 

The counter term for the dilaton depends on its conformal dimension and thus the 
quadratic term in V(ip), but we will not need its precise form here. For completeness 

11 



we carry out the holographic renormalisation of the background in full generality in ap- 
pendix |XJ Remember that we did not specify the details of this fall-off, only that it be 



at least constant or faster. Let us now expand the action (2.1) including the counter 



terms to quadratic order in the fluctuations. After integrating by parts the Euclidean 
renormalized on-shell action S os + Sct reduces to the finite boundary term 



W\&,4S>'' 



lim 

z-*0 



dojcPx 

2tt 



J2(n®)9tMz)4 

i 

e) 



*)(- 



u,z 



-E^4°(w^)4 0/ (-«^) 



Ak 2 z 



2~2 



g tx (uJ,z)g' tx (-UJ,z) - —g tx (u,z)g tx (-u,z) 



(3.12) 



where we have already Fourier transformed the time variable and used the fact that 
<&(z) — > 1 at the boundary. The condition that the solution be asymptotically AdS means 
that f(z) -> las well. For the configurations we are interested in there is no dependence 
on the field-theory spatial directions and so the integral j d 2 x = A will simply give us an 



overall factor of volume which explains why we included it in the definition (1.2). The 



fluctuations of the gauge fields and the metric have the near-boundary expansions 



Afi^z) 
9tx 



A^°(u)+zA^ l + . 



(0) 

9tx 



3 9tx "T" 



(3.13) 



Solving these equations with ingoing boundary conditions at the horizon will relate the 
subleading coefficients to the leading ones, defining the coefficient matrix C: 



4 } iM = E c ^ 



,w 



)4 i} '» 



C iP (u;)g tx (uj) 



9? x \u) 



C PP ^)gfX^) + Y, C ^ A ^ 



w 



(3.14) 



The matrix Cy is symmetric in its indices and time-reversal invariance implieaj that 
r - 2 9 2 r 

^Pi — -p-^iP- 

The generating functional will give rise to the desired correlation functions of the dual 
field theory by differentiating with respect to the sources. However, there is still a residual 
bulk diffeomorphism invariance as the condition g zx = only partially fixes the gauge. 
Since we started with a gauge invariant action, we must also have that under the trans- 







(3.15) 



formation (3.10) the function W is invariant: 

6 ies W[g° tx ,AV>°] 

6 Paraphrasing a standard argument [35j: Both P and JW are current operators, so they are odd under 
time reversal. Assuming that the vacuum is invariant under T we must have Gpj(ui) = rjjrjpGj'p^) = 



12 



This implies relations between some of the coefficients in (3.14), which have the interpre- 



tation of Ward identities. As noted below (3.10) we can also think of this transformation 



as resulting from translations of the dual field theory and so the resulting Ward identities 



are a consequence of this translation invariance. One finds that (3.15) implies that 

4^ 



C 



pp 



<e>, 



C 



Pi 



t 2 



(n w ) 



(3.16) 



Upon differentiating the on-shell action we obtain 

2tt 



A 



(P(wi)P(w 2 )) = ((e) + (p))S(u 1 + u 2 ) 



2n 



A 



(P{u x )J®M) = (n^)S(u 1 +u 2 ) 



2vr 
~A 



(j^M^M) 



~Y^ ^M + C JiM\ *( w i + <*>*) 



(3.17) 



Differentiating a generating functional fixes the retarded correlation functions up to pos- 
sible contact terms. We added a contact term proportional to the pressure (p) to the first 
line, which follows from translation invariance [36, 37J. All dynamical information about 
the conductivities is now encoded in the frequency dependent coefficients dj(u). With 
let us define the heat current 



j Q = p-J2^j (i) 



(3.18) 



(JqJq)- 



(3.19) 



Then the transport coefficients are 

Here a^ is the conductance matrix, describing both miralayer and mferlayer conductance. 
R is the thermal conductivity, which in our system is universal to both layers and finally 
a,i is the thermoelectric coefficient for each layer. We obtain the following simple relations 
between transport coefficients 



TR 



Tot; 



iyw«+ 



i ,3 

i 




2j> 



W( n W) 



ijj 



( n W) - ^3) a .. 



(3.20) 



The identities (3.20) generalise the relations of [20] to the double-layer case, albeit without 
considering magnetic fields. Clearly it would be interesting to extend our results to include 
them. Care must be taken when comparing to the discussion in J3S] as the B — )■ and 



to — y limits do not commute. We shall now analyze the transport coefficients (3.19) at 
low frequencies using the memory matrix approach, which we briefly review. 
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4 Low-frequency limit and momentum relaxation 

As mentioned several times already a translation-invariant system at finite charge density 
has a delta function divergence in its dc electrical conductivity. It is convenient to separate 
this divergence from the regular part, by writing the conductivity matrix as 

a ij (cj,T) = D«{T) f I + 7c5(u)^} + a^(u 7 T) . (4.1) 

We refer for D %: >(T) as the Drude weight, since it becomes the area under the Drude peak 
once translation invariance is broken. This divergence is associated with the conserved 
momentum. In the case at hand the divergent part of the conductivity is given by 

D«(T)= { l } } n J , (4.2) 

(e) + (p) 

so that each matrix element of cr 4 - 7 has a delta function divergence at zero frequency. 
The more precise reason for this divergence is that each current Jw has a non-vanishing 
overlap with the momentum operator, that is (J^'P) ^ 0. This implies that a finite- 
current configuration carries a finite momentum (P) = -^-(Jj). Since momentum is 
conserved this component of the current 'parallel' to P persists at arbitrarily late times 



and results in the delta-function divergences in (4.2). If translation is weakly broken, so 
that we can think of P as an approximately conserved quantity, the result is a Drude-peak 
in the conductivity with Drude weight given by D l: >(T), as we shall see in more detail 
below. 

Now observe that the Drude weight cancels for precisely the combination cr dra s intro- 



duced in (3.3). This is equivalent to the statement that (J-P) = 0, i.e. that the current 
mediating a dc has no overlap with the momentum operator. This implies that a current 
J_ does not have a component 'parallel' to the momentum and is thus able to relax on its 
own intrinsic timescale. The arguments presented here are somewhat heuristic (although 



the result (4.2) is precisely true). We now give a more rigorous derivation of these results 
using the memory matrix formalism, which will also allow us to make the stronger state- 
ment that <r dc is independent of the the mechanism of translation breaking, so long as it 
can be characterized by a coupling of the clean system to the leading irrelevant operator 
breaking translation. Before we proceed to the calculation itself, it is necessary to intro- 
duce a certain amount of formalism, more details of which can be found, for example in 
[38| [TO] . We will use the notation of the former reference for the remainder of this paper. 
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4.1 Some background on the memory matrix approach 



The memory function approach to transport, to the extent that we use it here, is a 
reformulation of the usual Kubo approach, suitable for considering late-time effects per- 
turbatively [39]. In this approach one develops a perturbative formulation for self-energies 
in correlation functions, in situations where a perturbation theory of the correlation func- 
tion itself is ill defined. We introduce, for sake of keeping the paper self contained, the 
essential notions here. 

We start by defining the inner product of two operators A and B at inverse temperature 
(3 = l/T 



A 



B 



) = /3 _1 / {A\0)B{i\))d\ 



(4.3) 



where ( • ) denotes an equilibrium average. We also have the time-dependent correlation 
function 



Cab® 



A 



-iLt 



B 



(4.4) 



where L is the Liouville operator, LO = i [H, 0\. This inner product satisfies the same 
hermiticity and completeness properties as the standard inner product between a bra and 
a ket. Denoting a general set of operators by O n , we note that the inner product defined 
here allows one to express their susceptibilities in a simple fashion: 



Xmn{T) = /3(O ro |0„) 



(4.5) 



The Laplace transform of (4.4), denoted Cab(w), is related to the retarded response 



function. More precisely we can use it to write the Kubo formula (1.2) in the form 



^( U )=fidj iJj {u). (4.6) 

Formal manipulations, similar to the ones involved in deriving the fluctuation-dissipation 



relations, establish that the Laplace transform of (4.4) can be expressed in the form of a 
dispersion relation 

Cij(uj) = iTXim [uXmn + iM^u)]' 1 \nj (4.7) 



with 

M mn {u)=P[0 l O,. 
where we defined the projector 



Q 



u-QLQ 



Q 



d t O n 



Q = l-Y\Om)px£n(p r , 



(4- 



(4.9) 
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So far we have done nothing more than to rewrite correlation functions and to repackage 



the Kubo function in a different form. If the sum in (4.9) runs over all operators in 



the theory, then the second term in (4.9) is a resolution of the identity and nothing is 
gained. However, as we shall see momentarily, the whole point of the formalism is that 
it allows us to focus on only a subset of operators, namely those operators which have 
a non- vanishing overlap with the (approximately) conserved quantities of the system. In 
the case at hand, these are all those operators that have a non-vanishing overlap with 
the momentum operator P. One often refers to these as the slow modes, as they relax 
on scales much larger than the typical microscopic scale of the theory. One then defines 
the projector Q by summing only slow modes, and the formalism becomes a very useful 



tool to study the late-time physics. The purpose of writing (4.7) is to repackage the 
relevant physical content of Cij{uS) in terms of a different function M mn (u). This memory 
matrix has a good perturbative expansion at late times, where a perturbative approach 



to Cijiuj) is impractical due to divergences of the form (4.2). The susceptibilities Xmn in 



(4.7) incorporate the information about the fast modes relevant to the late-time response. 



The form of the correlation function (4.7) is reminiscent of a dispersion relation and in 
particle physics language M(uS) would be a matrix of self energies as alluded to above. 
We will now bring the formalism introduced in this section to bear on the double-layer 
system which is the subject of our study. 

4.2 Double-layer memory matrix and DC limit 

Henceforth we assume that the effect of translation breaking is mediated by an irrel- 
evant operator which we add to the original Hamiltonian of the system, causing the 
non-conservation of the momentum via the Heisenberg equation of motion 

H = H + g5H => p = t [H 1 P}= ig[8H, P] . (4.10) 

This can be achieved either via umklapp scattering off a lattice, in which case 5H = 
G(kL), with lattice scale fc^, or it could be more generally via a disorder potential 5H = 
J dkV(t,k)0(k). In any case we assume that the coupling g is irrelevant, so that the 
system still flows to the same infra-red fixed point. 



The relevant slow modes, i.e. the modes going into a low-frequency expansion of (4.6) 



are the set {Ji, J2,P}- These then form the set of quantities {O n } entering the memory 



matrix (4.8) and projector (4.9). Thus considering the zero-frequency limit of (4.6) and 



(4.7) one arrives at the dc expression 



r ijf \ XpjWXpjU) 1 /a 11\ 

hma 3 (u) = -, (4.11 

w^o Xpp -L 
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where we should use our previous results Xpp — lim^o {PJ }■> Xpp — hm a; _ Sh o(-P-P) f° r 
the susceptibilitieq^J and the scattering rate T is determined by the PP matrix element 



of the memory function. Using (4.10) this boils down to the two-point function of the 



translation breaking operator, for example [TO 

r = ? 2 MM lim I 



L lim —lm.Goo{u-,kLj 



9=0 



(4.12) 



in the case of umklapp scattering. In the case of random disorder, the inverse rate r = 1/T 



is often referred to as the impurity timescale and the analog of (4.12) is derived by 



considering some averaging over disorder and thus involves an integral over momentum. 
The operator O would be the density operator J t for electric impurity scattering or 
scattering off an ionic lattice. 



Substituting the expressions (3.17) for the static correlation functions obtained above 



into this formula leads to a finite dc conductivity with Drude weight equal to the weight of 



the delta function in (4.2). Moreover, if the translation breaking operator is not present, 



that is if c? — > 0, then r — > and the divergence in (4.2) is recovered. 



4.3 DC limit of cr dra s 

Recall now that we noticed that the matrix element <T drag (o;) = o (u) has a finite dc 

limit even in the translationally invariant case, and moreover that its value is universal 
for holographic CFTs at finite density. We now consider the effect of translation breaking 



on this quantity. One can use the expression (4.7) for the this case, but since now the 



relevant current J_ has no overlap with the approximately conserved momentum there is 
no advantage in using the memory matrix. Its relaxation is entirely determined in terms 
of microscopic scales contained in Xmn(w). So in fact it is easier to go back to the original 



Kubo function (4.6). We find the answer 



lim a (u) 



lim — ImG R _ 

m0 oj 



(w, k = 0) 



9=0 



(4.13) 



to leading order with corrections parametrically suppressed by powers of the irrelevant 
coupling g. We recovered the same expression for <r dc that was obtained previously in the 
clean limit. This is perhaps not so surprising since we are explicitly assuming that the 
effect of the translation breaking is irrelevant in the sense of RG, so one would naively 
say that this must also imply that the effect on conductivities is negligible. But we wish 
to remind the reader that even such soft translation breaking has a dramatic effect on 

7 Comparing with [33 [TU] there is an order of limits issue: we should take to — > first and then B — > 0. 
I thank Chris Herzog for correspondence on this point. 
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charge transport in systems where the current operator overlaps with the almost conserved 



momentum. This is evidenced by the form of the low- frequency conductivity (4.11) for 
a typical element of the conductance matrix. It is because a dc is chosen such that it has 
no overlap with the conserved momentum, i.e. Xj~p = 0> that its relaxation is given in 



terms of the microscopic timescales of the CFT as computed in (4.13), and not in terms 



of the hydrodynamic timescales associated with conserved quantities of a typical element 



of the conductance matrix (4.11). 



It is important to emphasize that the arguments in this section have nothing to do 
with the holographic setup, relying only on the fact that there is no net momentum 
transfer involved in the transport process. In this sense it seems somewhat similar to the 
observation of Damle and Sachdev [133 that charge transport near quantum-critical points 
is finite at the particle- hole symmetric point. With this in mind one could also check 
that (the sum of) the Azlamasov-Larkin-type diagrams contributing to a perturbative 
treatment of drag (see FigjlJ) for the processes contributing to cr dc are insensitive do 
translation breaking in the sense above. 

The absence of a Drude peak in a drag is caused by the vanishing cross-susceptibility, 
Xpj- = 0. Let us briefly return to the case where only a Z 2 layer-exchange symmetry 
remains. Clearly P is even under this Z2, as is J + , which must be suitably defined as 
a sum of the two currents. However, J_, being a difference of the two currents, is odd, 
so that P and J_ cannot mix in correlation functions. This then ensures that the cross 
susceptibility Xp.J- vanisheaj Hence again we get a vanishing Drude weight and a finite 
dc limit of cx dras . 

We conclude this section by returning to the original definition of transconductance 



in section 1.1 Recall that this is related to induced charge flow the passive layer, when 
applying a current in the active layer. The boundary condition on this measurement is 
that the current flow in the passive layer is zero. It is a simple matter of algebra to use 
the general relations Jj = (f-'Ej together with the zero-current constraint, to show that 

E 2 =L 2 - a ^y 1 J 1 := P ^J 1 , (4.14) 

where the second equality defines the drag resistivity coefficient p drag . The Drude form 



(4.2) implies that the quantity in the round brackets is finite and well-defined in the 



zero-frequency limit, with the Drude weight canceling between the two factors. The zero- 



current condition together with the Ward identities, as distilled in D l ^{u) of Eq. (4.2) 



3 This Z2 is broken unless {rv- 1 ') = (n^), but one can let the Z2 swap the densities as well, i.e. define 



a weighted difference of currents as in (3.3 1. Nevertheless the quick argument here should be read with 



healthy skepticism unless the densities are equal. 



implies that the delta function contributions cancel. This cancellation is a close kin of 
the lack of Drude weight in the subtracted heat conductivity k, defined by measuring 
heat flow with a zero-current boundary condition. This was recently investigated in the 
context of a memory-matrix approach in [12] , where it was pointed out that k is similarly 
a universal quantity for CFTs. 

5 Discussion 

Let us conclude by evaluating the main points of this paper and their connection to 
other works. We reiterated the point that charge transport is a subtle physical process 
[SI El El EU HH [13 and therefore a careful analysis is needed when looking for universal 
values in the context of holography. Since charge transport usually involves momentum 
transport, and since momentum relaxation typically happens at hydrodynamical scales, 
the low-frequency behaviour of conductance (and thus that of resistance) is governed 



by hydrodynamical modes. We have seen that the quantity a dc defined in (3.3) involves 
charge transport without momentum transport, and as such does not have a Drude weight 
associated with it. Its dc limit is thus determined by the microscopic scales of the strongly- 
coupled field theory and we have seen that it is indeed universal. It would be interesting 
to analyse this quantity using a Boltzmann approach or another suitable technique where 
one can clearly see the interplay of hydrodynamic and microscopic scales involved. Note 
that spin conductivity (see e.g. [26^ 40J) is a quantity similar to our drag conductivity and 
one might hope that similar remarks about dc transport apply there. Going away from 
the dc value, we argued that the optical conductivity <7 drag (a;) is independent of frequency 
for certain strongly-coupled CFTs at finite density. Recently there emerged numerical 
[HJ S2J [H] evidence that strongly-coupled theories with a gravity dual can show robust 
scaling laws in the mid-IR region of the optical conductivity. Similar results have been 
obtained analytically by considering massive gravity in the bulk [43], albeit with scaling 
exponents that depend on unfixed parameters. In each of these cases the starting point 
was a translationally-invariant setup with a non-zero Drude weight, i.e. a conductivity 



of the form (4.2). It would clearly be interesting in this light to investigate the optical 
properties of our er drag (a;) as an example of a system with vanishing Drude weight at 
leading order, firstly in its own right, and also because it might give us some clues about 
physical processes in bad metals, which similarly have zero Drude weight [33]. 

The second half of this paper considered the effect of momentum relaxation on the 
interlayer conductance at the lowest frequencies within the formalism of Mori's mem- 
ory matrix. We pointed out that a generalization of the argument in [TO] implies that 
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each component of a n '(w) has a Drude weight determined entirely in terms of the Ward 
identities for translation invariance. This furnishes an alternative derivation of the hy- 



drodynamic form (4.2) at the lowest frequencies. Since the conductance matrix is the 
outer product of the charge vector (rii) with itself, there is always a linear combination of 
conductivities for which the Drude weight cancels and we showed that this combination 
corresponds precisely to the transport coefficient <7 drag . We also showed that scattering 
off a lattice as well as weak disorder has a negligible effect on this quantity, so that the 
clean value persists to the case without translation invariance. It is thus tempting to 
speculate that a dc is a candidate quantity to characterize intrinsic charge transport of a 
strongly-coupled critical theory. It seems then to be very desirable to investigate in more 
detail if certain physical systems, such as graphene heterostrucures, can serve as the ba- 
sis to fashion the requisite strongly-coupled layer systems to compare to the holographic 
computation. Of course we should not forget that the Coulomb interaction mediating the 
drag in real graphene (see for example jl5]) is in fact instantaneous, whereas our setup 
is fully relativistic. We nevertheless hope that this work more generally motivates the 
study of interlayer effects in strongly-coupled theories via holography, with the possibility 
of comparing to drag measurements. In this context it would be very interesting to study 



the transconductance in the presence of an applied magnetic field and include the Hall 

ij 
' ab 



components of cr* J 
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A Holographic renormalization 

Let us now be more specific about the holographic renormalization of the backgrounds 
we used implicitly in the main part of this paper. We gave a complete description of 
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the fluctuation part there and ignored the dilaton field ip as it does not play a role in 
the fluctuations. It does however play an indirect role via its backreaction and thus its 
influence on the thermodynamics and conserved charges. With this appendix we file the 
required paper work. Let's assume that the potential V(ip) has quadratic term m 2 . Then 
standard results mean that it is dual to an operator of conformal dimension 

m 2 -A^(A^-3)=0. (A.l) 

We require that the operator be normalizable, so that < A < 3, i.e. that the dual 
operator be relevant - or at most marginal. Its near-boundary behavior is then given by 
the expansion 

<p(z) = <^ A (1 + ■ ■ ■ ) + V2Z 3 ~ A (1 + ■ ■ • ) (A.2) 

with two integration constants (p\ and <p2- The naive on-shell action then diverges and 

we must add a local counterterm 

3- A 



qip 



J v^ 2 • (A-3) 



CT ~ 2kH 

If the conformal dimension is in a certain window [UJ there is an alternative quantization 
scheme, which amounts to adding a second counterterm and swapping the definition of 
source and expectation value, but we shall ignore this issue here. We can now compute the 



holographic stress tensor by varying the entire action (3.11 ) with respect to the boundary 
metric 7^. The result is 

k 2 (T^ v ) = lim UCp, - ICj^ - -7^ - -^-p 2 T|B ,J , (A.4) 

where K^ u is the extrinsic curvature of a constant z hypersurface (with outward pointing 
unit normal) and /C is its trace. One finds that 

K 2 (T%) = ^(A-3)^ 2 , (A.5) 

so that the theory is conformal in the UV if the product ^1^2 vanishes, or if the operator 
is exactly marginal, i.e. A = 3. In that case the trace of the stress tensor vanishes and 
the pressure is half the energy. We find it convenient to define 

lim -Too = (e) - ^~<Px<P2 (A.6) 



with 



and hence the pressure 



2£ 2 y /(1-V7) 

— -hm 



(e) = — lim M „ v " (A.7) 



(p> = ^ - -zzrvm ■ ( A - 8 ) 



(e|_3_-A 

2 2k 2 

Again, we see that in the conformal case the pressure is determined in terms of the energy 

by the zero-trace condition on the stress tensor. 
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